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In this study, we have considered the three dimensional mixed ferro-ferrimagnetic ternary alloy 
model of the type AB p Ci_ p where the A and X (X=B or C) ions are alternately connected and have 
different Ising spins S A =3/2, S s = l and S c =5/2, respectively. We have investigated the dependence 
of the critical and compensation temperatures of the model on concentration and interaction param- 
eters by using MC simulation method. We have shown that the behavior of the critical temperature 
and the existence of compensation points strongly depend on interaction and concentration param- 
eters. In particular, we have found that the critical temperature of the model is independent on 
concentration of different types of spins at a special interaction value and the model has one or two 
compensation temperature points in a certain range of values of the concentration of the different 
spins. 
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I. INTRODUCTION 



II. THE MODEL AND ITS SIMULATION 



Molecular-based magnetic materials have recently at- 
tracted considerable interest and study of the magnetic 
properties GLjlJl Jl 0, i, 0, & 1 EE El El EE 
El EE fl^.Tl7rfli. (ia |. A special class of the these 
materials, the so-called Prussian blue analogs, such as 
(Xj J Mn{i p ) 1 . B [Gr IW (CN)a].nH 2 (X"=Ni",Fe /J ) 0, 
5] and (Ni^Mn^Fe^ / )i. 5 [Cr /// (CN) 6 ].nH 2 @] which 
exhibit many unusual properties, for instance, occurrence 
of one [4] or even two jy] compensation points, mag- 
netic pole inversion [E 0] > the photoinduced magnetiza- 
tion effect [1, H , inverted magnetic hysteresis [l(| • These 
ternary alloys have ferromagnetic-ferrimagnetic proper- 
ties since they include mixed both ferromagnetic (J > 0) 
and antiferromagnetic (J < 0) superexchange interac- 
tions between the nearest-neighbor metal ions. The the- 
oretical investigations of these systems are difficult be- 
cause of their structural complexity. However, to ob- 
tain magnetic properties of the molecular-based mag- 
netic materials, up to now, these systems have been stud- 
ied by using effective-field theory [ll[ , mean field theory 
11411 and Monte carlo simulation (MC) methods 

in 

In this study we consider three dimensional ferro- 
ferrimagnetic ABpCi_p ternary alloy, consisting of three 
different Ising spins A=3/2, B=l, and C=5/2, which 
corresponds to the Prussian blue analog of the type 
(N^ / Mn{i p ) 1 . 5 [Cr /// (CN) 6 ].nH 2 j|. In this system, 
the coupling Cr-Ni is ferromagnetic and Mn-Cr is anti- 
ferromagnetic. Our aim, in this study, is to clarify the 
effects of the concentration and the interaction parame- 
ters on the magnetic behavior of the three dimensional 
ternary alloy model by using MC simulation method. 
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Three dimensional ferro-ferrimagnetic AB p Ci_ p Ising 
model consists of two interpenetrating cubic sublattices 
as seen in Fig. 1. It can be assumed that the A ions are 
located on the first cubic sublattice and the B and C ions 
are randomly distributed on the second cubic sublattice 
with the concentration p and 1 — p, respectively. Also, to 
construct a Hamiltonian for this system, the ion A can 
be represented by spin S A , and on the other hand, ions 
B and C can be represented by Ising spins S B and S c , 
respectively. If the interactions between nearest neigh- 
bors can be chosen such as A ions ferromagnetically in- 
teract with B, on the other hand, antiferromagnetically 
interact with C ions, thus, spins of the Prussian blue ana- 
log of the type (Ni^Mn{L p )i. 5 [Cr /// (CN) 6 ].nH 2 can be 
represented by this model where S A , S B and S c corre- 
spond to Cr, Ni and Mn, respectively. In this study we 
also consider next-nearest neighbor interactions between 
spins S A . 

The Hamiltonian of the considered system can be writ- 
ten in the form 

H = - J2 Sf[J AB Sfe 3 + J AC Sf{l - 

<nn> 

-j AA J2 s f s k (i) 

<nnn> 

where S A = ±3/2, ±1/2 for A, S B = ±1,0 for B and 
S c = ±5/2, ±3/2, ±1/2 for C, on the other hand, e 3 is a 
random variable which takes the value of unity if there 
is a spin X (S B or S c ) at the site j, if it not is zero. 
In Eq. (1), the first sum is over the nearest-neighbor and 
the second one is over the next-nearest neighbor spins. 
In this Hamiltonian the nearest neighbor interactions are 
chosen as J A b > and J A c < 0, and the next-nearest 
neighbor interactions are chosen as J AA > 0. 

In order to show the effects of the concentration p and 
the interaction parameters on the compensation and crit- 



ical temperature of the three dimensional ternary alloy 
model, we simulate the Hamiltonian given by Eq. (1). 
To simulate this model, we employed Metropolis Monte 
Carlo simulation algorithm [2(| to the L x L x L three- 
dimensional lattice with periodic boundary conditions for 
L = 10, 12, 16, 20, 24. One of the cubic sublattice is fully 
decorated with spin S A , and spins S B and S c are ran- 
domly distributed on the other cubic sublattice with the 
concentration p or 1 — p, respectively. All initial spin 
states in the L x L x L three-dimensional lattice are ran- 
domly assigned. Configurations are generated by making 
single-spin-flip attempts, which were accepted or rejected 
according to the Metropolis algorithm. To calculate the 
averages, data, over 20 different spin configuration, is ob- 
tained by using 50000 Monte Carlo steps per site after 
discarding 10000 steps. 

The sublattice average magnetizations per site are ob- 
tained by 




where Nb denotes the number of B ions Nb = pL 3 /2, 
whilst Nc represents the number of C ions Nc = (1 — 
p)L 3 /2 on the same cubic lattice. Total magnetization 
per site is given by 

M = i (M A + M B + M c ) • (3) 

III. RESULTS AND DISCUSSION 

In this section, we have given the simulation results of 
the ternary alloy model AB p Ci_ p and we have also dis- 
cussed the dependence of the critical and compensation 
temperature on the concentration and other interaction 
parameters in the Hamiltonian. Simulation results have 
been obtained for the system with lattice size L = 10, 
12, 16, 20 and 24, however, here, we have only presented 
the results of the model with lattice size L = 20. We also 
note that the critical temperature of the system for the 
different interaction rates and concentrations have been 
obtained by using of the method of the finite-size scaling 

m - ra- 

In a recent study [151 ] it was reported that two dimen- 
sional ternary alloy model does not show a compensation 
temperature point when there is no next-nearest neigh- 
bor interactions term in the Hamiltonian i.e., Jaa = 0. 
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FIG. 1: The crystallographic structure of prussian blue analog 
with two interpenetrating cubic lattices. 
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FIG. 2: Dependence of the critical temperature on interaction 
ratio 7? in the three dimensional ternary alloy AB p Ci- p for 
different values of p when Jaa = 0.0. 



However, our simulations show that the system has a 
compensation point for all R (we set R = \Jac\/ Jab) 
values in interval of 0.1 < R < 2.642 at p = when 
Jaa — 0. This point will be considered below. Now, 
in order to compare with the previous results [TEI, [HI], 
in Figs. 2 and 3 we discuss the dependence of the critical 
temperature of the three dimensional ternary alloy model 
ABpCi_p on interaction rate R and concentration p for 
next-nearest neighbor interactions i.e., Jaa = 0. 

In Fig. 2, the critical temperature of the three dimen- 
sional ternary alloy model has been plotted as a func- 
tion of R for various values of p when Jaa — 0. It can 
be seen from Fig. 2 that the critical temperature of the 
system has a linear dependence on the interaction ratio 
R and there is a critical behavior at a special R value. 
When R c = R = 0.513, the critical temperature of the 
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FIG. 3: Dependence of the critical temperature on the con- 
centration p in the three dimensional ternary alloy AB p Ci_ p 
for several values of interaction ratio 7? when Jaa = 0.0. The 
lines show a part of the second-order transitions separating 
the ferrimagnetic and paramagnetic phases. 



system has a fixed value of T c = 5.47 for all p values. 
At i? c , the critical temperature of the system does not 
change with concentration p. This means that neither 
the spin-I ions nor spin-5/2 ions substitution to system 
change the critical temperature of the system at R c . This 
critical behavior has been reported in theoretical and ex- 
perimental studies [HI, [H, [3 ■ The value of the R c for 
ternary alloy AB p Ci_ p whose spins consist of S A = 3/2, 
S B = 1 and S c = 5/2 has been obtained as R c = 0.4781 
in the study based on mean field approximation [l8| and 
as R c = 0.49 in the Monte Carlo simulation of two dimen- 
sional system [l5j ]. Furthermore, the experimental mea- 
surements indicate that there are Prussian blue analogs 
at the R = 0.45 have a T c which is almost independent 
of p [HI . Fig. 2 also reveals that concentration p plays 
an important role for the ternary alloy model AB p Ci_ p 
since it determine the kinds of the spins and interactions 
in the system. For example, when p = 1 and p = 0, 
the system AB p Ci_ p fully reduces to the ferromagnetic 
mixed spin-3/2 and spin-1 and ferrimagnetic mixed spin- 
3/2 and spin-5/2 Ising system, respectively. As seen in 
Fig. 2, although T c of the system is independent of p at 
R c , however, the total magnetization of the system may 
considerably change owing to relatively small variation of 
the concentration p. Indeed, for different values of p, the 
dependence of critical temperature of the system on the 
interaction ratio R is very different above and below of 
R c . This behavior can be explained by the change of the 
concentration p in the system. On the other hand, it can 
be detected from Fig. 2 that when R < R c , the critical 
temperature of the mixed spin-3/2 and spin-5/2 system 



FIG. 4: Magnetization of the three dimensional ternary alloy 
ABpCi_p vs temperature for different values of p (Jaa = 7.5 
and R = 1.0). 

is smaller than mixed spin-3/2 and spin-1 system. On 
the contrary, when R > R c , the critical temperature of 
the mixed spin-3/2 and spin-5/2 system has the highest 
value. On the T c lines, the critical temperature of the 
mixed spin-3/2 and spin-1 Ising system is equal to that 
of the mixed spin-3/2 and spin-5/2 Ising one. 

In Fig. 3, the dependence of the critical temperature of 
the three dimensional AB p Ci_ p system on the concen- 
tration p has been shown for several values of R when 
Jaa — 0. The lines represent part of the second-order 
phase transition separating the ferrimagnetic and para- 
magnetic. Fig. 3 provides the argument that the concen- 
tration p determines the magnetic features of the system 
mentioned above. Indeed, Fig. 3 clearly shows that the 
critical temperature of the system is changed by the con- 
centration p for fixed values of R. As seen from this 
figure that, when R < R c , the critical temperature of the 
system linearly increases with increasing of p, whereas, 
when R > R c , the critical temperature of the system lin- 
early decreases with increasing of p for fixed values of R. 
However, when the values of R close up R c , the criti- 
cal temperature of the system more slowly, but linearly, 
change with increasing p, and at the critical R c value, the 
critical temperature of the system denoted by triangle- 
line in Fig. 3 is independent of the concentration p. On 
the other hand, Fig. 3 also shows that the interaction rate 
R plays an important role on the critical temperature of 
the three dimensional AB p Ci_ p system. Finally we state 
that the critical temperature of the model are consistent 
with previous result (l8| . 

In this study we recognize that the three dimensional 
ternary alloy model AB p Ci_ p has one compensation be- 
havior for Jaa — 0, however, for Jaa 7^ it has one 
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FIG. 5: Magnetization of the three dimensional ternary alloy 
ABpCi_p vs temperature for different values of p (Jaa = 7.5 
and R= 2.642). 



FIG. 6: Dependence of the compensation temperature T C om P 
on the interaction parameters in Hamiltonian of ternary alloy 
AB p Ci_p for p = 0.25. 



or multi compensation points, when other conditions are 
satisfied. However, the appearance of the compensation 
temperature is strongly affected by the interaction and 
concentration parameters. Indeed we see in the present 
study that the model has not a compensation point for 
all values of p and R. The dependence of the compen- 
sation temperature behavior on concentration and other 
interaction parameters has been discussed below. For 
discussion, although the system has been simulated in 
the intervals of 0.0 < p < 1.0 and 0.1 < R < 2.642, the 
value of Jaa used in the present study is chosen based on 
previous theoretical study [l5[ • The results of simulation 
for R — 1.0 and R = 2.642 are respectively represented 
in Figs. 4 and 5 for chosen parameters. 

One compensation point has been found in the in- 
tervals of 0.0 < p < 0.3 and 0.1 < R < 2.642 when 
Jaa — 7.5. However, it is seen that the system has not 
compensation behavior for the same values of parame- 
ters when p > 0.3. For R — 1.0 and several values of 
p, the compensation behavior of the system can be seen 
from Fig. 4. On the other hand, as seen from Fig. 5, the 
considered system has a multi compensation behavior at 
R = 2.642 and p = 0.3 for Jaa = 7.5 while it has one 
compensation point for 0.2 < p < 0.3. Furthermore, our 
simulation data introduce that the system shows com- 
pensation behavior at p = for all values of R, when 
Jaa = 0. In addition, in the case Jaa — 0, the com- 
pensation point has been found for R = 0.25 at p = 0.2, 
0.25; for R = 0.75 at p = 0.3; for R = 1.25 at p = 0.3; 
for R = 2.0 at p = 0.1, 0.3, 0.4; for R = 2.642 at p = 0.1, 
0.2, 0.3, 0.4. Whereas, it has been reported in previous 
study that there is no compensation point for Jaa = 



in two dimensional model [151 ] . 

The effect of the interaction parameters on the com- 
pensation behavior of the three dimensional ternary 
model is also discussed in Fig. 6. This figure shows de- 
pendence of the compensation temperature T comp of the 
model on interaction parameters in Hamiltonian only for 
a fixed value of the concentration parameter p (p = 0.25). 
In this figure, square-line represents the behavior of the 
compensation point vs Jaa for fixed values of Jab = 5, 
J Ac = ~5 and p = 0.25, circle-line indicates the behav- 
ior of the compensation point vs Jab for fixed values of 
Jaa = 7.5 and J ac = — 5 and p = 0.25, and on the other 
hand, the behavior of the compensation point vs J ac 
is plotted for fixed values of Jaa = 7.5, Jab — 5 and 
p = 0.25 with triangle-line. As seen from Fig. 6 that for 
fixed p, Jab = 5 and Jac = — 5, the compensation tem- 
perature decreases slowly as the strength of the Jaa in- 
creases. Similarly for fixed p, Jaa — 7.5 and Jac = — 5j 
the compensation temperature decreases slowly with in- 
creasing of Jab- However, for fixed p, Jaa = 7.5 and 
Jab — 5, the compensation temperature dramatically 
increases as \Jac\ increases. These results indicate that 
the compensation temperature has a strong dependence 
on the parameter Jac whereas its dependence on Jaa 
and Jab is relatively weak. The characteristic behav- 
ior of the dependence of the compensation temperature 
on the parameters of present model consistent with the 
results of two dimensional model [15| . 
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IV. CONCLUSION 

In this study, we have considered the three dimensional 
ternary model AB p Ci_ p whose spins consist of S A = 3/2, 
S B = 1 and S c = 5/2. We have investigated the depen- 
dence of the critical and compensation temperature be- 
havior of the considered model on concentration and in- 
teractions by using MC simulation method. We have ob- 
served that the behavior of the critical temperature and 
the existence of compensation points strongly depend on 
interaction and concentration parameters. Particularly, 
we have found that the critical temperature of the model 



is independent on concentration of different types of spins 
at a critical R c value and the model has one or two com- 
pensation temperature points in a certain range of values 
of the concentration of the different spins. We concluded 
that magnetic properties of the system AB p Ci_ p can be 
controlled by changing the relative concentration of the 
different species of ions. As a result, we would like to 
stress that these theoretical results can be very useful 
for designing molecular magnets in experimental studies 
since the existence of compensation in the ternary alloy 
ABpCi_p that can be setup by adjusting the proportion 
of compounds. 
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